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PEEFACE TO THE FIEST EDITION. 



The work now oflFered to the University is strictly an 
Elementary Treatise. No attempt has been made to go 
into all the varied details, of methods and examples, 
which present themselves in the wide field of Partial 
Differential Equations, considered purely as an Algebraical 
subject. 

I have endeavoured, however, to omit no important 
consideration affecting the Principles of those Equations. 
And I trust that the methods of solution here explained, 
and the instances exhibited, will be found suflScient for 
application to nearly all those important problems of 
Physical Science, which require for their complete in- 
vestigation the aid of Partial Differential Equations. 

G. B. AIRY. 



Royal Observatobt, Greenwich, 
1866, August 15. 



PEEFACE TO THE SECOND EDITION. 



Seyebal verbal alterations are made in this Edition; two 
small paragraphs are added; and some sentences are 
introduced, referring to works in which examples of the 
application of the Theory of Partial Differential Equations 
may be found. But nothing is changed in the plan of 
the work, and no alterations are made in the numbering 
of the Ai-ticles. 

G. B. AIRY. 



BoTAL Obsebvatoby, Gbeenwich, 
1873, July 15. 
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ON PARTIAL DIFFEKENTIAL 

EQUATIONS. 



FRELIMINABY NOTICE ON INTEQBATI02C. 

1. In all that follows, we Bhall suppose that it is 
always possible to effect simple integration ; inasmuch as 
any difficulties of integration, connected with the solutions 
of Partial Differential Equations, do not affect the principle 
of those solutions. Thus, for instance, we shall not hesi- 
tate to represent an unknown function of x by x' (a?), 
(the second differential-coefficient of xi^))> ^^ *^® assump* 
tion that, whatever be the form of x{x), we can in some 
way find the function ^ (^) of which it is the second diffier- 
ential-coefficient. 

CHARACTERISTICS OF THE SOLUTIONS OF SIMPLE 
DIFFERENTIAL EQUATIONS. 

2. Before entering on the subject of Partial Differ- 
ential Equations, it may be convenient to consider some of 
the characteristics of the solutions of Simple Differential, 
Equations. 

3. To begin with Simple Differential Equations of 
the first order. Suppose, for facility of geometrical illus- 

tration, we consider the equation y -^ = a, or the equation 

B. E. B 



2 ON PARTIAL DIFFERENTIAL EQUATIONS. 

dx , 

y ^ = h,oi. which the former, translated geometrically, in- 
dicates that the subnormal of a plane-curve (to be found) 
is constant, and the. latter indicates that the subtangent is 
constant. The algebraical solutions are easily found : in 
each, there is a constant, which does not occur in the 
original diflferential equation, and \b not defined by it ; a 
constant of that class described (perhaps improperly) by 
the term "arbitrary," but which really means "not yet 
determined, but enabling us by proper determination of 
its value 60 to fix the value of x corresponding to a given 
value of y that we can adjust the solution to some specific 
condition/' 

The reader is requested to observe that instead of the 
term " arbitrary constant," we shall always use the term 
" undetermined constant." 

4. If we treat the geometrical translation of the dif- 
ferential equation by a geometrical process, always draw- 
ing a normal or a tangent (as the case may be) so as 
to make the subnormal or subtangent constant, then 
drawing by means of it a small portion of the curve, then 
repeating the process, &c., we may produce a polygon 
which will approach to the strict solution, with smaller 
errors (by taking the sides small enough) than any small 
quantity that can be assigned. In each case, however, a 
starting-point is necessary. • 

5. In both ways (the algebraical and the geometrical) 
of treating the problems these conditions manifestly 
hold :— 



c . 



CHAEACrrEillSTICS OF SIMPLE BlFFERKln^IAL EQtTATIONS. S 

The curve, in each case, is one definite curve. 

The curve, in each case, is ^ continuous curve, ex- 
pressed by the same equation through - its T^hole 
extent. 

Even if there be isolated points or curves, still the 
same one equation defines the whole, 

fl 

It is necessary to introduce one undetermined quantity^ 
enabling i^s to adjust the curve to a specific con- 
dition given by special considerations: that unde- 
termined quantity is however a simple constant. 

6. Let us now consider a Simple Dififerential Equation 
of the second order : such, for instance, as is given by this 
problem, " To find the curve in which the radius of curva- 
ture is a function, given in form, of the ordinate." Here 
the algebraic solution gives a formula requiring two unde- 
termined quantities, still simple constants : the equivalent 
geometrical treatment shews that we require two elements 
(as, for instance, the value of x and the inclination of the 
tangent, for some one value of y). But all the conditions 
hold which are mentioned in Article 5 : the only difference 
being that, instead of one undetermined quantity, a simple 
constant, there must now be two undetermined quantities^ 
simple constants. 

7. The conditions which will be found to hold in the 
•solutions of Partial Differential Equations differ very re- 
markably from these. 



B2 



4 ON PARTIAL NFFEKENTUL EQUATIONS. 

« 
INTRODUCTION Or TWO OR MORE INDEPENDENT 

VARIABLES. 

8. It is convenient, in Simple Differential Equations, 
to contemplate one quantity z as expressed algebraically 
in terms of another quantity x (whether it be actually so 
expressed, or not) : and to consider its differential-coeffi- 
cients 

djs d /dz\ d^ o 

dx^ dx\dx) da?* *' 
as being formed by taking the limiting values of fractions, 
in each of which the denominator is an increment in the 
value of X, and the numerator is the corresponding incre- 

that X is the " independent variable.'* Geometrically it is 
illustrated by supposing that we consider an ordinate z of 
a plane curve, as also the inclination-tangent of the 
curve's-tangent, the change of that inclination-tangent for 
a small change in the value of x, &c., to be expressed in 
terms of the abscissa x, 

9. There is no difficulty in conceiving z to depend on 
two quantities x and y, combined in any way and with any 
constants under any functional formula; and in considering 
that we may at pleasure change the value of x without 
changing that of y, or may change the value of y without 
changing that of x, or may change both simultaneously. 
If we do not vary y, y is pro tempore a constant : and by 

dz d^z 
varying x alone, we may form -7-, 1-^, &c., just as in the 



TWO VARIABLES IK SOLID GEOMETRY. 5 

functions which we have illustrated by reference to a plane 
curve. If we do not vary x, then x is pro tempote a conf- 

stant, and by varying y alone, we may form -j- , -jjf , &c. 

But if we vary both x and y, then we have the two series 
of differential-coefficients which we have just set down, 

and also -^ j- or -j j- (which, as is known in the ex- 

dx ^ ay ay ^ ax ^ 

pansions obtained by the Differential Calculus, have the 

same value), with other coefficients of analogous form, 

which do not appear in the succession of coefficients, such 

as those of Article 8. Here we consider x and y as " two 

independent variables," which are strictly independent 

of each other : and £? as a function of both. 

10. Algebraical Geometry of three dimensions assists 
greatly in illustrating these algebraical conceptions. If a?, 
y, IS be three rectangular co-ordinates, then the expression 
of z by means of aj and y determines the numerical value 
of the height of an ordinate z that must be erected over 
the point on the plane ajy, which is defined by any nume« 
rical values of x and y, in order that the elevation of the 
summit of that ordinate z may represent the value of our 
algebraical function. If we do this for an infinite number 
of values of x and y, we determine an infinite nximber of 
ordinate-summits, aQ which he in one curved-surface. If, 
supposing this curved-surface formed, we then contem- 
plate the values of z with x invariable and y variable, we 
shall include all these values which are in a plane parallel 
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to the plane yz, and at the given distance x from the 
origin of co-ordinates. This is the same as forming a sec- 
tion of the curved-surface by a plane parallel to the plane 
yz. and considering the curved-intersection as an ordinary 

plane-curve^ Thus we can ofctain -j- , -j-^ , &c. as in Ar- 

dy ajf , 

tide 8 : but it must be remembered that, in general, x is 
lurking in their expressions as a constant ; and therefore, 
in the complete diflFerentiation (with respect to a;) of any 
formula in which they may occur, each must be considered 
as a function both of x and of y. Similarly, the values of 
z with y invariable and x variable will be formed by a sec- 
tion of the curved-surface, made by a plane parallel to xz : 

the differential-coefficients -y- , -^ , &c. can be formed ; 

and a similar caution with regard to their dependence both 
on X and on y must be borne in mind. 

11, The coefficients -j- , -7--, &c. are, as in Article 8, 

ax flwr 

to be considered as defining the inclination- tangent of the 

curve's-tangent, the change of that inclination-tangent for 

a small change in the value of x, &c., &c., in the curve 

formed by the intersection of the plane parallel to xz 

with the curved-surface- In like manner, the coefficients 

T' > T^ > &<5. define similar elements in the curve formed 
ay ay* 

by the intersection of the plane parallel to yz with the 

cPz 

curved-surface. But the coefficient -5 — 7- requires further 

ax. ay ^ 



TIME IS OFTEN ONE VABIABLE. 7 

explanation. Its real expression is -i- [ -r-J • Now -j- de- 
fines the inclination-tangent of the curve's-tangent in. the 
sectional plane parallel to yz : and therefore its differential- 
coeflScient with respect to x defines the rate at which that 
inclination-tangent, in the plane parallel to yz, is changed 
by shifting the sectional-plane in the direction of a?. 
Analogous explanations apply to analogous succeeding 
terms, 

12. Algebraically, there is no diflSculty in conceiving 
i3 as a function of any number of independent variables, 
as % V, w, X, y, and in treating the various diflferential co- 
efficients, which may rise to any degree of complexity; 
But we cannot extend any further the elucidation derived 
fi:om Solid Geometry. 

13. Although, in strictness, no difference is really 
made by any physical peculiarity in the nature of the 
quantities represented by the independent variables, pro- 
Tided that their magnitudes can be defined by numbers 
and can therefore be treated by algebraical process'; yet it 
may be well to mention that, in some of the most important 
applications of Partial Differentials, one of the iiidependent 
Variables is the expression for time. Thus, suppose that 
we are considering the nature of the disturbance which' 
each particle of air is undergoing in a musical pipe : we 
want to know the entire movements of the particle whose 
original ordinate was X; and this implies that, in the 
general formula for the disturbance, x (the general symbol of 



8 ON PARTIAL DIFFERENTIAL EQUATIONS. 

X) is to have the special value for that particle, and ia 
therefore to be made constant for that particle, but that no 
limit is to be put on the variations of L We also want to 
know the state of disturbance of all the particles at a 
Certain tinie T, and this implies that t is then to be made 
constant, but that (r is to admit of all possible variation. 
And these can only be comprehended in a general formula 
which admits of variations both of x and of t Similarly 
for the motion of the particles of waves of water, &c., &c. 

14 In Simple Diflferential Equations, the data of 
the problem, whether mechanical or geometrical, usually 

lead to an equation between os, is, -r- , -y-^, &c., from 

which we desire to obtain a general expression for z in 

terms of x. Similarly, in Partial Differential Equations 

the data of the problem usually lead to an equation be- 

dz dz d^z d^z d^z « #• i • i 

tweenx,y,z.;^, ^. ^, ^_^ ,_, &c. from which 

we desire to obtain a general expression for z in terms of 
w and y. It is the object of the present treatise to shew 
how this may be done, and how the meaning of the solu-, 
tions may be explained, in some of the simpler cases : and 
we now proceed with the Solution of Partial Differential, 
Equations of the First Order. 
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TREATMENT OF THE SIMPLEST PARTIAL DIFFERENTIAL 
EQUATION OF THE FIRST ORDER. 

15. The simplest Partial Differential Equation of the 
first order is 

dz dz dz dz ^ 

dan ay ax ay 

[Here, and in all subsequent operations, we shall put 
V for flw? + y.] 

Evidently this equation is satisfied by making z^u^ 

For T- = o* -7- = Ij aJid therefore -3 — a -7- = 0. But 
dx ay ax ay 

it is also satisfied by making z = ^(ii)> where <^ expresses 

a function whose form is totally unlimited. For then 

dz ±t f \ du ,f, s 

and therefore 

g-a| = f(«).{a-a} = 0. 

16. We have thus found a solution of considerable 
generality in form ; but there is no evidence that the 
generality is complete. In order to obtisiin complete and 
certain generality, we shall use a process of Change of 
Independent Variable, guiding ourselves in the selection of 
a new Independent Variable by the indications derived 
from the last imperfect process. 
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17. Instead of considering js as a function of x and t/, 
let us consider 2; as a function of x and u. [Such a sup- 
position is certainly competent to represent z, because, if 
we had z expressed by x and y, we have merely to sub- 
stitute for y its value u^ax, and z then becomes expressed 
by X and u.'] Having z tlien a function of x and u, where 
u is itself a fdnction of x and y, we proceed to express the 
original equatipn in tenns of the new differential-coeffi- 
cients. We shall put ~-^ for the entire differential-co- 

efficient in regard to every way in which x could appear 

dz 
when X and y were used, and -j- for the differential-coeffi- 
cient when X and u are used: the original equation is 
therefore to be written 

ax ay 

18. 2; is a function of x and u, where u depends on 
X and y. 

Therefore ■' 

d{z) _ dz dz du ^ dz dz 

dx dx du' dx dx du * 
and 

dz ^ dz du ^ dz ^ 

dy du' dy du * 

Substituting these in the equation which is to be solved, 

d(z) dz dz 
ax • uy ^ ax ....... 
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That is to say, when z is caressed in terms ofx cmd u, 

dz 

— =0. 
dx 

The integration of this simple quantity involves the 
most iD(iportant considerations of the whole theory. 

19. By the ordinary rules of integration, z = G, where 
C is something which does not contain x. But what may 
it contain ? It may contain every thing whatever except x. 
It may contain any constants whatever. It may contain u : 
for the diflferentiation with respect to x will not touch u, 
which is another Independent Variable and is a constant 
with respect to x. It may contain any function of u. It 
may contain these included in the form 

" any function of constants and u" 

But, in ordinaiy algebraical language, this would be 
suflSciently described as " any function of u" provided that 
we bear in mind that u may be combined with any con- 
stants, and may even actually disappear, so . as to leave 
constants only. This being understood, the solution of the 
equation is 

z = ^(tt), 

or 2J = <]>{(ix + 7/), 

where <f> expresses a function of any form whatever which 
considerations not yet presented to us may induce us to 
adopt. 

The steps of this process leave no opening foi: further 
generality in the solution : which is therefore complete. 
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20. Instead of assuming z to be expressed by x and w, 
with elimination of y, we might have assumed z to be 
expressed by y and u^ with elimination of x^ Then the 
• process is 

d(z) ^ dz dz du ^ dz dz - 
dy dy du' dy dy du ' 

dz dz du dz 



dx du' dx du * 



and 



^ dz ^ d{z) _^ dz ^ 

dx dy "~ rfy' 

from which the same conclusion follows, namely 

= ^(w) = (f>{ax + y). 

21. We may even change the Independent Variables 
80 as to eliminate both x and y ; and, as this process will 
be important hereafter, we shall shew its application here^ 
Let u=:ax + y, v = ex+Jy, where e is not equal to af; 
and let z be supposed to be expressed in terms of u and v. 



Then 



^^ — ^ ^ d^f dv _^ dz j_^^ 

dx du' dx dv' dx "^ du dv ' 

dz __ dz du dz dv dz -t ,^^ /, 

dy "" dw * ^ dv' dy ^ du dv • ^ 
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and the equation 

dz dz ^ 

dx dy 
becomes 

/ A\ dz ■ yv dz A 

whence, in the same manner as before, 

GEOMETRICAL INTERPRETATION OP THE SOLUTION. 

22. In order to shew the geometrical meaning of this, by 

reference to such a curved-surface as has been considered in 

Articles 10 and 11, we may observe that if hx and Sy are 

exceedingly small increments of afand y, then the increment 

A dz ^ dz ^ , . , . dz dz y 

ot = -7- oa; +-y- ov ; which, smce -r- = « -7- * becomes 
dx dy ^ dx dy 

dz 

-T- {pZx + hy) ; and, if we always make Sy = — a8a?, the 

increment of z will be = ; that is, if we always pass from 
one point on the plane xy to another point on that plane by 
a short line making with the axis of x the angle whose 
trigonometrical-tangent is a (which process continually 
repeated will produce a finite straight line in the same 
direction), the value of z will be invariable. Therefore, the 
curve-surface will consist of a series of parallel lines, in 
each of which the values of z are the same so long as 
ax •\- y\A the same, but which is subject to no other con- 
dition. It will therefore be such as is represented in figure 
\\ a cylindrical surface (in the widest sejise of the word), 
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whose axis is parallel to xy and makes with x the angle 
whose trigonometrical-tangent ia «, but whose transversal 
section, or section by either of the planes xz, yz, is abso- 
lutely undefined by the equation. Those sections may be 
adapted to specific physical data, and are therefore of the 
nature of the quantities which we call undetermined. But, 
as the only condition at present established is that the lines 
(of which the surface consists) shall be parallel, the sections 
need not to consist of curves defined each by a single 
equation in its whole extent. A section may consist of a 
bit of a circle with a bit of a cissoid and then a bit of 
a logarithmic spiral; and these may be joined at any 
angles ; or the section may be utterly undefinable by alge- 
bra : and after it has gone on to some extent, there may be 
an absolute interruption and then it may begin again, &c. 



CHARACTERISTICS OF THE SOLUTIONS OF PARTIAL 
DIFFERENTIAL EQUATIONS. 

23. Thus it appears that these conditions hold in the 
solution ; 

It is not certain that one symbolic expression will 
represent the solution, or that one curve defined by 
one equation will represent the section of the curve- 
surface. 

It is not certain that the section is continuous, either as 
to actual connexion of parts, or as to gradual change 
of direction, or in any other circumstance. 
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^iThere may be isolated parts of the curve, of totally 
distinct character, 

I 

An undetermined quantity must be introduced : but this 
undetermined quantity will not usually be a simple 
constant, but will be a function of x and y, and (in 
the most general case) a discontinuous function, such 
as corresponds to the discontinuity of curve mentioned 
in Article 22. 

On comparing these conditions with those in Article 5, 
the reader will see the characteristic differences between 
the solutions of Simple and of Partial Differential Equa^ 
tions. 

INTERPRETATION OF THE SOLUTION WHEN ONE 
INDEPENDENT VARIABLE IS TIME. 

24. There is another instance of the solution in 
Article 19, &c. or (more frequently) of a solution of the 
same form in the treatment of Partial .Differential Equa- 
tions of the Second Order, which merits explanation. Let 
one of the Independent Variables be time : put t for y, 
and suppose the ftinction of t + ax to be -^ {i (< + ax)} or 
'\lt{bt + ex), c being = ba. Then z = y^{bt + ex). Let any 
one value of bt + cx he JB: then if t be increased by cG 
(where C is really arbitrary), and if x be increased by 
^ bC, the value o{ bt + ex is still = J9, and therefore z 
remains the same. That is, at a time later by cG, we shall 
have the same value of z, provided that we take an abscissa 
smaller by bC. This is just the same as if the ordinate z 
was slid backwards upon the abscissa x with a velocity 
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represented by -^, or - , The same applies to every other 

(ordinate z. And thus it appears that the whole curve 
defined by the summits of z may be conceived as changing 

•L 

its shape by sliding backwards with a definite velocity - . 

c 

This is the characteristic of a wave, in which (as in waves 
of water, to take the most familiar instance) the form 
travels while the particles do not travel with the form. 
Whenever therefore we meet with such a term as ^ {ax + y) 
we may at once take for granted that one of its applica- 
tions may be to the motion of a wave. 

If the value of z had been yjr(bt'' ex), we should have 
found in like manner that the wave is travelling forwards. 

These remarks do not in any degree interfere with 
those of Article 23. 



OTHER PARTIAL DIFFERENTIAL EQUATIONS OF THE 

FIRST ORDER. 

uiS uZ 

25. We proceed with the equation ^ — a-r- =^a{x, y), 

where the form of the function a is given. We shall solve 
it by the method of Change of Independent Variables. 
Let u = ax + y, v = ex +fy, where e is not equal to of, 
(which, by giving different values to e and /, includes the 
three changes in Articles 19, 20, 21).-^'Trhe value of x in 

terms of u and v ia'^ — , that of y is — ^ . Treat- 

af-e' ^ af-e 
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ing the equation as in Article 21, we find 

whence, integrating on the same principles as before, 

where, in the integration with respect to t?, w is to be con- 
sidered constant; and, after the integration, ax'\-y and 
ex -^fy are to be substituted for u and v. 

26. On trying this process upon any ordinary function 
(as a function in integral powers of x and y, or a circular 
or exponential function) of which the integrations are easy, 
it will be found that e and /disappear from the result. It 
may be well to shew this in a simple instance. 

Let one term of a {x, y) be a;',y'. 

The corresponding term of the integral is 
-1 



' J V 



I 

Integrate by parts, beginning with the first factor; it 
gives 

(p-m.t.z-.r'^' O-'''"-'-""^'"')' 

P. K. C 
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and, by repeating the process, we shall obtain the result in 
5 -f 1 terms, in powers of x and y, and without e and /. 

27. If we had begun the integration-by-parts with 
integration of the second factor, we should have obtained a 
different series of terms of powers of x and y ; yet the 
difference would have been such as to produce no difference 
in the final result,, as applicable to any physical investiga- 
tion in which the undetermined function is to be adapted 
to given physical circumstances. A simple instance will 
shew this. 

Let a{x,y) = X, j^. Then we obtain the two following 
series for the integral : 

First, 
z = ^{ax+y) +^^-2{6aVy' + 4aVy + aV}. 

Second, 

1 

z = ^jr{ax + y)+j^{-4axy'-y% 

The second value, it is easily seen, may be put under 
the form 

t {a^ +» + j^, {6aV/ + 4a Vy + a V} 
- j^ {/ + 4air/ + 6a'a;y + 4a Vy + a V} 
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^ yJr {ax + y) - j^, [ax +y)* 
+ j|^ 16aV/ + 4aVy + aV]. 

Now — :j-^-3 (a« + yY is a function of {ax + y) ; and, when 

we proceed to determine the form of the function y^ (ax +y), 
which must be adopted, in order, with the other terms of 
the expression, to satisfy any condition prescribed by 
physical considerations, we must at the same time take 

into account — ^^^ (ax+y)* as a function of (ax+y). We 

•may therefore at once unite it with y^ (ax+y). Call their 
sum X (^^ + y)* Then the second solution is 

« = X (aa; + y) + j2^, (6aVy' + 4a Vy + a^x*) ; 

precisely the same as the first, except that the symbol x 
stands in the place of <f>. But in both cases the form of 
the function <^ or ;^ will be determined by the consideration 
of satisfying the same physical condition. And therefore, 
whether we write ^ or ;)^ we shall afterwards arrive at the 
same expressions; and therefore the two solutions, fir^t 
and second, are in use identical. 

28. In many cases, when the form of solution has 
been obtained by the process of Article 25, the details will 
be obtained most readily by assuming the form of solution 
with indeterminate coefficients. 

c2 
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29. It is not intended here to go into any discussion 
of the more complicated forms of Partial Differential Equa- 
tions of the first order. The following however may be 
mentioned as flowing immediately from the solutions which 
we have found. 

and the equation is immediately reduced to the form 
already treated. 

dss dss 

If the first side is € (a?) x -^ — ax ^(y) x -j- , 

and the equation is reduced to the same form. 

dz^ ^ dz f f > dz^ ^ dz 
dx " dx dy ^ dy* 

and the equation becomes 

^i-a— * « 0. 
dx dy 

In reference to physical investigations, the theory of 
Partial Differential Equations of the first order is princi- 
pally valuable as introductory to that of the second order. 
In this view, it is hoped that the instances here given are 
sufficient. 
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TREATMENT OF THE SIMPLEST PARTIAL DIFFERENTIAL 
EQUATION OF THE SECOND ORDER. 

30. The simplest Partial DifiFerential Equation of the 
seoond order is 

dx * dy 

Integrating with respect to a?, and remarking that, by 
virtue of the reasoning in Article 19, there must be added 
to the integral an undetermined function of jr. 

Integrating this with respect to y, and remarking that, 
for the same reasons, an undetermined function of x must 
be added, 

GEOMETRICAL INTERPRETATION OF THE SOLUTION* 

31. This equation and its solution admit of easy- 
geometrical illustration. Referring to Article 11, it will be 

dz 
seen that -v- is the trigonometrical-tangent of the angle 

between y and the curve-tangent, in the curve formed by 
the intersection of a curved-surface with a plane parallel to 

zy. And the equation -^ — -r- ^ 0, or -i- ( -r- ) = 0, denotes 

dz 
^ that, when x is changed, j- undergoes no change ; or that 
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the inclination in the plane ty will be the same for a 
point more or less advanced in the direction of a; as for 
the point under consideration. The curve-tangent may be 
absolutely higher or lower in the direction of z, but. it will 
preserve the same inclination to y. As this applies to 
every point of the intersection-curve, it follows that the 
intersection-curve more or less advanced in the direction 
of X will have the same form as that at the point under 
consideration ; or that all sections by planes parallel to zy 
give the same curve, though perhaps at different elevations. 
And from this it follows (considering that, when two similar 
curves are separated in. the direction of x, the slope from 
every point of one to the corresponding point of the other, 
in the direction of a;, must be the same) that all sections by 
planes parallel to zx give the same curve^ but not necessa- 
rily the same as those given by planes parallel io zy. Then 
the elevation of any point z may be considered as composed 
of these two parts ; first, the elevation of the corresponding 
point of that zy curve whose plane passes through the 
origin, which elevation (since the form of the curve is the 
same for all values of x) may be called ^ {y) ; secondly, the 
elevation of the point z above the point last considered » 
which elevation (since the points are connected by a curve 
whose form is independent of y) may be called -^ (a?). 
The whole elevation of the point z will therefore be 

32. These functions may be discontinuous, for th^ 
reasons stated in Article 22. 
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33. In figure 2 is a representation of one of the solids 
in which ^ {y) and ^ (x) are both continuous functions. 
It is scarcely necessary to remark that the forms of such 
soUds may be infinitely varied ; thus (taking the vertex as 

origin of coordinates), if js = — hT-, or both sections be 

par^ibolas with different parameters, the curve-surface is a 

paraboloid with elliptic base; if a = — = — , the curve* 

surface is the paraboloid of revolution (from which, figure 

a? V* 

2 was drawn); if 5? = v, the curve-surface is like a 

ah 

limited portion of the interior of an annulus, or like a 

mountain-pass, &c &c. &c. 

In figure 3 is a representation of a curve-surface in 
which one function is continuous (the sections being similar 
parabolas), the other is discontinuous (each section being 
two sides of a triangle, the triangle being the same 
throughout). 

. Ill figure 4 is a representation of a surface in which 
both functions are discontinuous (each section being two 
sides of a triangle, but the triangles in the plane xz being 
different from those in the plane yz). It forms, in fact, a 
pyramid with trapezpidal base ; two vertices of the trape- 
zium having the same value of x, and two having the 
same value of y, and the vertex of the pyramid being above 
the intersection of the diameters of the trapezium. 
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34. The solution of the equation -7 — -y- = a (a?, y) 
obviously is, 



z 



= ^(y) + ^W+j .j .a(a?,y). 



TREATMENT OF ANOTHER PARTIAL DIFFERENTIAL EQUA* 
TION OF THE SECOND ORDER. FIRST, BY CHANGE OF 
INDEPENDENT VARIABLES. 

35. To solve the equation 3;;^"^'"^ = ^ (^> y)' 

[This equation is the most important of all, especially in 
reference to those physical theories in which wave- 
transmission, of sound, or of light, or of water, &c., is 
explained by mechanical theory.] 

The best method of solving this equation is by the 
Change of Independent Variables. 

Let t£ = aa: + y, t;==aaj— y, 

/ , . , . u-Vv u — vX 

^which givea?= -2^,y ^-y j. 

Consider as a function of x and y because it is a 
function of v, and v. 

Then 

dz _^ dz du dz dv ^ dz dz 

dx du dx dv dx ^ du dv ' ' 

d^z _ d (dz\ _ d (dz\ du d (dz\ d'd 
do^ dx \jdIx) "" du \dx) dx dv \dx) dx 
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d fdz dz \ d fdx . dz \ 

=5 jT j-x« +:r-x«P<« +j-(t7-x^ •^-ry<a]xa 
au \au av J av \au av J 

da' d«* . dv rft;' ' 

dz ^ dz du dz dv ^ dz -i .dz ^ ^ 
dy ~^ du dy dv dy du dv ' 

d^z _ d (dz\ d (dz\ du d (dz\ dv 

dy ~ dy \dy) du \dy) dy dv \dy) dy 

d fdz ^ ,dz ^^ -, , d fdz ^ dz -\ ^ 
di* Vdi^ dv J dv \du dv J 



r«. 



di*' du . dy dw" ' 



Tx tt <2f .d 5? A 9 d z 

Hence -r^ — » -r"2 =* 4a-; r* 

oa? d^ dw » dv 

And the original equation, divided by 4a', becomes 
d'^ 



du . dt; 



^ 1 fu A-v ti — 1?\ 



Whence, by Article 34, 



or 5? 



. / V . / \ 1 f f /tt + V «— 1?\ 
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where, after the integration, for u and v are to be put ax-\-y 
and ax^y. 

35* If the, right-hand term of the given equation 
were a simple function of a;, multiplied only by constants 
and the first power of y, as if 

g-«'|^ = (6 + ci^)xy3"(^). 

let z^ = z — {b + cy)x fi (x): 

then -^ = ^^(b + cv)xff'(x) —^ = — • 
and the equation becomes 



d\ 



--a'^/ = 0; 



dx^ dy^ 

of which the solution (by the formula above, when a = 0) 
is z^ '= ^ {ax—y) + '^{ax +y) ; 

or z ^ <f> {ax — y)+'yjr {ax -{-y) + {h + cy) xfi ix), 

35**. When a = 0, the solution « = ^ (w) + ^/r (v) may 
be represented geometrically in the same manner as 
z = <f>{x)+'^{y) in Article 33, with this difference only, 
that u and v are ordinates on the plane a;y, which are not 
at right angles, except in the case when a = 1. 

36. When one of the independent variables is time, 
it will be seen by the same reasoning as in Article 24 
that the two undetermined functions represent two waves 
travelling in opposite directions. 
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SECOND TREATMENT OF THE SAME EQUATION, BY SEPARA- 
TION OF THE SYMBOLS OF OPERATION FROM THOSE 
OF QUANTITY. 

37. A second method of solving the equation 

is founded upon the " separation of the symbols of opera- 
tion from those of quantity." This theory is, in fact, 
merely a convenient form for exhibiting the indubitable 
results of legitimate algebra ; but it sometimes serves to 
suggest new methods of treating equations, which, when 
verified, are useful. 

38. The equation 

dz dz . m^ -m-^ 

may be written, almost without departure from ordinary 
notation, 

d d 



{i*4>'(^-^-- 



the connexion of the left-hand bracket with z which follows 
it being, however, not by multiplication but by diflferen- 
tiation. And, if 

d{X, Y) cZ(X, Y) . . 



28 ON PARTIAL DIFFERENTIAL EQUATIONSi 

we may, on the same principle, write it 



(£-«|)(^'^)=«(^'^)' 



and therefore 



still on the understanding that the connexions on the left 
hand are not by multiplication but by differentiation. 



(; 



39. Now, if we treated all the left-hand symbols as 
signs not of operation but of quantity, their product would 
be 

or, pursuing the fanciful idea still further, 

and this suggests the idea that, if we follow up the true 
differential operations instead of the fanciful algebraical 
multiplications, we shall arrive at that result in its true 
differential meaning. And so, in fact, it proves. For, by 
actual differentiation, 

d fdz^, dz\ ^ ^*« , <?*g 
dx \dx dy) " da^ dx .dy^ 

^ /^ dz\ _ d^z %d^z ^ 

dy \dx dy) "" dx . dy^ dy*' 
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andj adding the two lines, 

( d ^ d\ fdz .dz\ J ^ _ i^Pz 
\dx ay) \jix dy) dx* dy*' 

Here the analogy with the algebraical operation has led to 
a real gain of convenience, by shewing that we can here 
l3reak up the operation for a Partial Differential Equatioi^ 
of the second order into two operations for the first order. 

40. Thus, to solve the equation 

d z ^d z . ^ 

we have first to solve the equation 

from which (^, Y) will be found ; and then to solve th$ 
equation 

from which z will be found. 

Both operations are effected by the process of Article 25; 
the most convenient assumption however being 

u^ax + y, v^ax'-y. 

We shall return hereafter to the " separation of symbola/' 



I 
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OTHER PARTIAL DIFFERENTIAL EQUATIONS OF THE 

SECOND ORDER. 

41. To solve the equation 
or 

It does not appear possible to solve this equation gene- 
rally, except by the use of imaginary symbols. If we take 
the second form of the equation (as written above), and 
apply the solution of Articles 35, &c. 

If, as is usual, we suppose ^ and -^ to represent real 
functions, it will be impossible to destroy the imaginaiy 
terms in this expression except by making ^ the same as 
<^, and thus giving up part of the generality of the solution. 
But the introduction of two undetermined functions may be 
thus obtained. Let 

where ;^ and q> are real functions. 
Then 
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Expanding the functions by Taylor's Theorem, it will 
be seen that the imaginary terms are destroyed, and there 
are two undetermined functions. 

It will be seen hereafter that, in some practical appli- 
cations of this equation under limiting conditions, the 
solution may be obtained more easily by specific process 
than by using the general solution. 

42. The following equation, which occurred to the 
author in the investigation of the movement of spherical 
waves of air (gravity being neglected), deserves attention 
for the nature of the substitutions made and the form of 
the result. 

The equation given by the physical considerations is, 

di? ' dx\x dx) ' 

The last term may be made simpler by assuming 

dx \xj ' 
or y ^ x\ z. 

J X 

r ■ 

Dififerentiating this twice with respect to a:, 



d^y « . dz 

and 1 ^ = ?f + ^. 

X dar X ax 
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whence the second side of the equation 



dx \x (Lx?) ' 



J2^ 

For the first side, or -r^ , we have 

av 

d'z cP d fv\ d fP 



de 



^ £_ jd (y\ ^ d_ ^ (y\ . 

d^ • dx ' W dx ' df \xj ' 



which, since ^ can be aflfected by ^^ only because y is a 
function of t, becomes 

dvKau'dfJ' 
and the original equation is changed to 

dx \x ' dt*J ' dx\x' da?) 



Integrating with respect to x, 

1 dV_o* ^ 
x' d^ X ' da? 



X" (<). 



Solving this by tlie method of Article 35*, 
whence 
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42* In the Author's Treatise on Sound, Article 43 
and several following articles, will be found a discussion of 
equations of the form 

where m has different integral values. 

43. We believe that scarcely any other equations than 
these, occurring in physical investigations, have been solved 
in a finite form. Several equations have been solved in 
the unsatisfactory form of infinite series, of which the con- 
vergence is not always assured ; but it is not the object of 
the present Treatise to enter on them. 

44. The following equations which have presented 
themselves in the author's investigations, (probably among 
many others occurring in physical inquiries), have not 
been solved. 

In investigating the vertical motion of a horizontal 
wave of air, on the supposition that the elasticity is propor- 
tional to the density, and that gravity is constant, we 
arrive at an equation of this form : 

d^X dX ,oCrX ^ 
djtr dx aif 
D. E. D 
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which we have pot succeeded in integrating generally in a 
finite forau 

In investigating the radial motion of a spherical wave 
of air whose center is the center of the earth, supposing the 
elasticity to be as the vl^ power of the density (where n 
must be greater than a certain quantity which is> 1, in 
order to make the whole mass of atmosphere finite), the 
resulting equation takes this form : 

( r )dr\ar rj r \ar rj 

This equation may be much simplified, but in any 
form we have not rendered it integrable. 



FURTHER CONSIDERATION OF THE SEPARATION OF THE 
SYMBOLS OF OPERATION FROM THOSE OF QUANTITY. 

45. In Article 39 attention was called to the method 
of solving the equation then under consideration by the 
''Separation of the symbols of operation from those of 
quantity." This principle may be applied to an extensive 
system of equations, which are however, at present, matters, 
rather of curiosity than of physical value. It may be 
sufficient here to indicate one class. 
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If the equation 
can be expressed as 

then the process of Article 40 can be applied with so little 
difficulty that it appears unnecessary for us to delay 
further on it. In the successive changes of Independent 
Variable, it will be found convenient to take the factors in 
successive pairs ; thus, for the effect of thq first pair ot, 
fiietors, where 

(s-|)(c-'|)«^=«('.A 

assume ax-^y^u, bx-^y^v, 

and proceed as in Article 25. 

46. There is one exceptional case> however, well 
brought to notice by this method of treatment, which 
merits further attention ; namely, that in which two factors 
Of b, are equal To take the simplest case, suppose 



ux axdy dy^ ' 



d2 
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or 



{A^a—^(—-a—]z = 
\dx dy)\dx dy) 



fliere we may assume ax'\-y = «, ex -{-fy s= v, e and / 
"being any constants whatever whose proportion is not a : 1; 
the value being admissible for either. 

Then 

dx du dv ' 

d (X Y) ^ d{X,Y) ^^ d(X, Y) , 
dtf du dv 

and, putting (X, Y) for (^- « ^) «» 

the equation (^ - a^j (X, Y) = 0. 

' / r\ d (X^ Y) ^ 

becomes (^ -«/)-' — ^ — = ^* 

whence {X, F) = ^ (u)^ 

And, putting for (X, Y) its value. 



(^-«|)^ = *(^^^ 



or, by the same substitution which was made in regard to 
the differential coefficients of (X, Y), 
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integrating this, and remarking that <f> (u) possesses the 
properties of a constant as regards integration with respect 
to V, 

or 

Observing that a constant factor may be supposed to be 
included in the form of an undetermined function, we may 
express z in either of the forms, 

y>^^i {y +^) •H'fL (y +«») ; 

all which, as applicable to the satisfaction of any assigned 
conditions, will be found to be identical, on employing the 
reasoning of Article 27. 

47. In pursuing this subject, the signs / and d may be 
considered as reciprocal, or / = d"*. In some cases, equa- 
tions may be varied by carrying symbols, such as 

dx dy' 

which hold the position of multiplier on one side of the 
equation, to the position of divispr oh the other side, qr 



S8 ON PAftTlAL DIFFEBENTIAL EQI^ATIONS. 

as multiplier in the form 



\dx dyj 



dy. 

And in this shape they may be substituted in abstruse 
and general theorems. For instance^ to take a case which 
(in this theory) is almost elementary, we know that 



(0(«'")-"(£-^i"e:-*^) 






(^ + **)"«• 



If n =s — 1, thia becomes, reversing the sides of the equation. 






But the equation 



dz dz _ 
dx dy"~ ' 



d . d 
or 



fd ^ d\ 



may, in accordance ivith these principles, be put in the 
form 



fd , dy* 
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and the second side agrees with the first side above, 
provided that ^t- he put for a and treated as a constant. 



Then we shall have 



2? = € dy 



J s 



a solution to which, in some cases, a real meaning may be 
given. 

48. This principle, as a purely algebraical and 
symbolical process, possesses very great power, and leads to 
very remarkable results. But the reader cannot fail to 
observe that it carries with it no evidence whatever for 
the validity of results (such as is conveyed by the opera- 
tions of quantitative algebra, or by the steps, properly 
pursued, of the diflFerential calculus), for which it must 
rely on subsequent verifications. As aiding the appli- 
cation of Partial Differential Equations to physical 
investigations it possesses little value. The further 
examination of it would therefore be out of place in 
this Treatise. The student who desires to follow it up 
will find much information in Boole's Treatise on Differ^ 
ential Equations, Gregory's Examples^ and similar works. 
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TREATMENT OF PARTIAL DIFFERENTIAL EQUATIONS WHEN 
THE SOLUTION IS LIMITED BY PREARRANGED CON- 
DITIONS. 

49. We have seen in Article 44 that there are some 
Partial DiflFerential Equations of which the solutions 
cannot be obtained by known methods, so long as the 
utmost generality is required in the solutions. But in 
some cases, by attaching limiting conditions to the solu- 
tions, the same equations may be integrated with ease; as 
will be seen in the following instances. 

50. In considering the tidal disturbance of water in an 
equatoreal canal round the earth, retarded by friction pro- 
portional to the velocity (see the Eacyclopoedia Metropoli- 
tana. Article Tides and Waves), we arrive at the following 
equation {x the original ordinate of a particle measured 
along the canal, X the horizontal disturbance of the par- 
ticle) : 

-W " ^«'^(^«-^) -f-dt +" dZ' 

where H, f, a*, are constants depending on the moon'3 
attraction, the coefficient of friction, and the depth of the 
sea. This equation cannot be solved generally. But, 
remarking that the only part of the solution which has any 
interest for us is that which follows the same law of pe- 
riodicity as the lunar motions, we may assume 

X = -4 sin {it —mx) + -Bcos [it — Tna?), 
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where A and B are constants to be determined. On 
substituting this, the values of ^ and B and the expression 
for X are found without difficulty. 

51. Before adverting to a specific solution of the 

equation -^ +a* -t-^ = 0, we will make the following 

remark on the general solution. When an exponential 
form is adopted for the functions in the general solution, 
Article 41, the general form may be retained conveniently: 
but in other cases it is usually necessaiy to expand the 
functions. Performing this process according to the ordi- 
nary algebraical rules, the general solution becomes 

-» W y +^ (y) f7273 " 

which, as will be found on trial, satisfies the equation. 
A solution by series, however, can rarely be considered as 
quite sufficient. 

52. In considering the motion of ordinary small waves 
in a sea of uniform depth {x horizontal in the direction in 
which the wave is going, y vertical measured upwards from 
the bottom), we find (see Tides and Waves) 

dy^ dx^ 

The cumbrous form of the general solution of thi^ 
equation would make it almost useless to attempt to apply 
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it to tlie special case. But the nature of the case permits 
us to assume 

X = -Bcos {nt-^mx) +8 sin {nt —wiar), 

where -B and 8 are functions of y. Substituting, we have 

whence 

and, after some reductions peculiar to the problem, 
X ^ A (€'"•' H-e"*") cos (nt-mx-B). 



FINAL DETERMINATION OF THE FUNCTIONS WHICH ARE 
UNDETERMINED IN THE GENERAL SOLUTION. 

53. In nearly the whole of these •solutions, our re- 
sulting expression for the quantity which it is our object 
to find is accolnpanied with undetermined functions. It is 
now an object of great importance to determine those func- 
tions. Although in reality the process for doing this may 
be stated in brief rules; yet it will be better understood 
from an exhibition of its application in one or two actual 
instances than in any other way. 

54. Example 1. In the consideration of the tidal 
motions of water without friction in an equatoreal channel 
round the earthy X being the horizontal displacement at 
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anytime f of a particle whose ordinate measured round the 
channel from a fixed origin was x, the equation is 

^ = jff sm (tt^mx) +c'-^ ; 

where it is double the moon^s hour-angle from the origin, 
and mx is double the longitude-arigle of the point x from 
the same origin; and where c' is a constant depending 
on the depth of the water. This equation can be com- 
pletely solved : its result is 

X - >. a .« sin (i<-maj)+^ (cb-^x) +^[ct-x) ; 

which we shall write 

Csin {it -moo) +<f> {d +x) +i^ (ct — a:)» 

Each of the three terms represents a wave, but they 
are waves of different characters. The first wave corre- 
spends in its velocity with the lunar force Hsin (it—mx) 
which produces it : for this the author has introduced the 
term /orcecZ wave; it depends entirely on the lunar force, 
and would not exist if there were no lunar force. The 
second wave is travelling in the direction opposite to 
the moon's apparent diurnal motion, and the third wave in 
the same direction as the moon^s apparent diurnal motion; 
for these the author has introduced the term free waves; 
their velocity is independent of the moon's motion ; the 
two waves, or either of them, may or may not exist ; they 
are entirely independent of the moon's action, and the 
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<luestion of their existence is generally independent of 
the moon's action and of the forced wave. This is the state 
of things when there is no material limit to the length of 
the canal. 

55. But suppose that the canal is limited by the ob- 
stacle of land at its two ends (as in the Mediterranean Sea^ 
all minor gulfs, &c., being put out of consideration). Then 
the forced wave alone is not sufficient for a solution. For, 
the existence of the terminal obstacles (whose ordinates we 
will call a and b) requires that, at those obstacles, the hori- 
zontal displacement of the water be nothing at all times ; 
and therefore when x = a or x = b, X must == whatever 
be the value of U 

Now this cannot hold with the forced wave alone ; for, 
the values which the first term assumes for a; = a and x = b, 
namely, (7 sin {it —ma) and (7 sin {it —mb), "Will not vanish 
for all values of t We must therefore have recourse to th6 
two undetermined functions; or, in other words, we must 
now necessarily introduce the two free waves as an indis- 
pensable part of the solution: and we must determine their 
elements so that the two conditions, of X always = when 
x = a and when a? = J, shall be satisfied. 

; 66. The term expressing the forced wave is a simple 
periodical term having it in its argument. Therefore, in 
order to destroy this at all times for definite values of x, 
the free waves must be expressed by simple periodical 
terms haying it in their arguments. Therefore the two 
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undetermined terms ^ {ct +a?) and -^ (c^ — a?) must have 
the forms 

A^ sin |- {ct +a?) 'j- +B^ cob j- (ct +x) I , 



and 



^.sin 1^ (c^-a:) J +5,cos|- (ct-x) j j 

and the entire expression for X will have the form 

C Bin (it -^msc) 

» • 

J" A^ sin (i^ -f - «r) +JB, cos (i^ +- x) 

c c 

+ -4, sin (i< — x) 4-5, cos (i^ — x) ; 

c c 

6r, expanding the trigonometrical terms, and putting e for 

I tX %0R I 

sin it X \ (7 cos mx -f c cos — + A sin — ^ 
I c c) 

+ cosi< X j — (7sm ma; +/sm — l-^'cos —V . 

Making this = when x — a and when aj = J, and^ 

«/* «^ 

writing a for — and /3 for — , we have the four equations, 

c c 



» (7cos ma +^ cos a +A sin a, 
- — £7 sin ma +/sin a +g cos a, 



I 
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= Ccosmb i-^coafi + hsmfi, 
=a — Csia m& +/sin fi +g cos /S. 

Solving these equations, 

e siri (/8 — a) = C (— cos ma . sin )S +cos mi . sin a), 
y sin 08 — a) = C (— sin ma . cos )8 +6in mb . cos a), 
g sin (J3 —a) = C ( sin ma . sin /8 —sin mi . sin a), 
A sin O— a) = (7(cosma,cos^-cos7ni.cosa). 

Substituting, and restoring the original notation, 

X = v-;j — ^ siJi U^ — mx) — . X 

(ow — i") sin j- (6 — a) >• 

sin (t^— ma). sin- -{i —a?)}- +sin(iit-m6).sin ]- (a? —a)!- L 

Thus the terms, undetermined in the general solution, 
are now completeily determined so as to satisfy the special 
conditions of the problem. The entire expression for X, it 
will be seen, satisfies the original partial differential equa- 
tion, and also makes X ~ at all times when X » a and 
fl? — i. 

57. Example 2. In the ordinary problem of vibrating 
musical strings, where no force is supposed to act after the 
string has been put in motion, x being measured longi-* 
tudinaJly along the string and z being the small transversal 
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displacement of any point (see Sound, Article 73), the 
equation for z is 

where a' = gL, L being that length of the same string 

"whose weight will represent the tension of the string : and 

its solution is 

z =s ^ (at +x) +'^ {at — aj). 

Here the solution is expressed entirely by the two 
undetermined functions. If the initial circumstances of 
disturbance (namely, the displacement and the velocity of 
every point) were given, those functions could be deter- 
mined, as will be shewn below. But even without an 
absolute knowledge of the forms of the functions, some 
very important properties of the solution may be ascer- 
tained thus. 

58. Let the ordinate of the near end of the string be 
0, and that of the further end I. Then, whatever be the 
value of t, z isO when a? = or when x = I. That is, 

<f> {at) 4--^ {at) « 0, 

^ {at +1) +1^ {at -Q = 0. 

The first equation shews that yfr must always = — ^, 
Thus one of the undetermined functions is already elimi- 
nated (the other will be eliminated from consideration of 
the initial circumstances). 

Substituting in the second equation, 

^ {at +1) -0 {at -I) « 0. 

Now at may, with changes of time, receive any value 
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whatever. And if at-^l = q (which makes at-^-l^q +2Q, 
q may have any value whatever. A^d thus, q having any 
value whatever, we always have the equation 

that is^ the form of the function ^ must be such that, on 
increasing or diminishing the quantity under the bracket 
by 21 f the same value of the function will be reproduced. 
Now in either term of the expressioti for the displacement 
of the particle at x, namely 

<f> {at +x) —<f> {at —a;), 

the increase by 2/ of the quantity under the bracket may 

21 
be eflfected by increasing at by 21, ot by increasing ^ by — . 

a 

21 
Therefore, every time that we increase f by -— , we find 

a 

the particle a: with the same displacement as before. And 

this applies to every particle of the stritig. Consequently 

the string performs a complete vibration, returning again 

. 2t 
to the same position, in the time — . 

It is true that, as will shortly be seen, under certain 
circumstances the string may perform two or more com- 
plete vibrations in the same time, but it always returns to 

21 
the same state after the interval — - ; and this is the only 

a *" 

interval of which it can be asserted that a return to the 
same state necessarily occurs in that time. 
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The reader may exercise his ingenuity in proving that 
the motion of the string consists in the simultaneous 
movements of two waves in opposite directions, each of 
which travels backwards and forwards from end to end 
with velocity a, and is reflected at each end. 

59. Putting, for the moment, Q for the quantity 
affected by the functional symbol <^, it appears that <^ {ff) 
is a finite periodical function, going through its changes 
and returning to the same value when 6 is changed by 11. 
A trigonometric function of sines and cosines of an angle 
goes through all its changes and returns to the same value 
when the angle is increased by 27r. Hence it appears that 
<^ (d) is similar in its general character to a trigonometric 

function of sines and cosines of — ^y— , or , Q. Any 

finite periodical magnitude may be represented with any 
assigned degree of approximation by a series of integral 

powers and combinations of the sine and cosine of — ^; and 

these powers may be converted into sines and cosines of 
multiple arcs. Thus the function ^ (5) may be represented 

by . 

A^ sin -^ Q +-4j sin -y- ^ +&c. + -4„sin -r- ^ + &c. 

+5jC0S^^ +jB,co!S -^^ + &c. +J5„cos — ^ + &c.; 
which may be conveniently expressed by 

D. E. £ 
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and thus the value of z will be 

2 I J. sin '^Z^i^l +S |5„ cos ^^E^+^X 

_2 1 A sin ^" ^"/ -^) } -2 15. cos !!!:i^)} , 



or 



Z I 2^^ cos -J- sm -y a; 1 



— S (2J3^sm — , SI 



WTT 



sm -. x]\ 



dz 
and the value of -,, will be 

at 



S^f2mra . . nirat . rnr \ 



-2 



/2r?7ra ^ nirat , nir \ 

( —J Bn cos — y sm -y iC J . 



60. Suppose now that the initial circumstances of 
displacement {Z) and velocity {Z') are given for eveiy 
point of the string, and that from these we desire to deter- 
mine all the subsequent motions. We must now make 
^ = 0, and we have 

Z = 2 [2^«siny a;j, Z' = — 2 f -y-njB„sin-y arj . 
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And OUT immediate object is, knowing the values of 
Z and Z for every value of a?, to determine the values of 
^j, A^ &c., B^y B^, &c. This will be done by pieans of 

the following general formula. If we multiply sin -y- a? 

by sin -y x, where m is any integer different from n, and 
integrate from a? » to a? = Z, we have for product, 

1 (m— w)w 1 (m+n)7r 
— cos ^ — ^—^ — X — TT cos -^ i— ^ — X ; 

for general integral, 

I , (m — n) TT Z . (m +w) ir 

— sm -^ f-^ — X — ,s-7 ^^ — sm ^ r-^— x; 

2 (m — n) TT f 2 (m +w) ir I 

and for definite integral, 0. 

But if m «s w, we have for product 

1 1 2n7r 



2"2'''''"T'^' 



for general integral 



X I . 2w7r 

2 -in^ir '"^ -r ^' 



and for definite integral, ^ • 
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Hence, if we put S for the integral with respect to x 
taken from a; = to a? = Z, 



SlZsin—j xj = A^l,&c. 

8 iZ' sin -y- a?) = — ttaw B^ ; 

8 iZ' sin —j- x\ = —man B^, &c. 



61. In a given instance, in which the initial displace- 
ments and velocities are given not by formula but by 
numerical values, the following process will apply. Sup- 

pose we limit the multiples of y 0, where sines and cosines 

are to be used, to the first six, namely. 



'^/) ^'^a^'^a ^a ^'^a ^'^/j 
jO, j^>iff> x^ T^ T^ 



which will suffice in almost any conceivable case. Divide 
the length I into 120 equal parts, and suppose that Z and 
Z' are known numerically for the middle of each of these 
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nir 



parts. For the first of these middles, -j- x = 45' ; for the 

second it is 3 x 45' ; for the third it is 5 x 45' ; and so on. 
And for each of these parts, the integral is, in fact, taken 

through the extent j^^ . Hence, the formula above gives 

the following numerical rule : — 

Multiply the successive values oi Zhy sin 45', sin 3 x 45', 
sin 5x45', &c.: the sum of products = 120 A^. 

Multiply the successive values of Z by sin 2x45', 
sin 6x45', sin 10x45', &c.: the sum of products = 120^,. 

Multiply the successive values of Z by sin 3x45', 
sin 9x45', sin 15x45', &c.: the sum of products = 120 A^, 

And so to A^. 

And a similar process (with proper changes) determines 
i9j, 5j, &c., from the given values of Z\ 

Then substituting these in the formulae at the end of 
Article 59, the circumstances of all subsequent motion are 
completely obtained. 



62. If A^, or JSj, or both, are found to have real 
values, and A^ A^ &c., B^, B^ &c., have no real values, 
then the vibration is that of a simple line of sines, de- 



54 ON PARTIAL DIFFERENTIAL EQUATIONS. 

pending, for time, on the angle —j- . If only A^ and JS^, 

one or both, have real values, then the vibration is that of 

27r 
two lines of sines ; separated at the point where -j-a? = tt, 

or where a? = ^ , thus forming a permanent node in the 
middle of the length ; and depending for time on the angle 

— ^ — , and thus vibrating in half the time of the former 

vibration. Similarly if A^ and jB, only have real values, the 
string is divided into three equal parts by two nodes, and 
the vibrations are made in one-third of the time of the 
fundamental vibrations. 

The connexion, of these results, with the results of the 
experience of the senses as to the musical tones produced, is 
one of the most important points in the Acoustical Theory 
of Music. 

4 

63. The two instances which we have given will 
probably suffice, better than any rules, to shew the use 
that is to be made of the undetermined functions in the 
solution of a Partial Differential Equation. For important 
applications, in the Theories of the simple Echo and of 
Resonance, see Soimd, Articles 41 and 81. 
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CONSIDERATIONS ON THE NECESSITY OF CONNEXION 
BETWEEN THE NUMBER OF UNDETERBONED CONSTANTS 
AND THE ORDER OF THE EQUATION (IN SIMPLE DIF- 
FERENTIAL EQUATIONS), AND BETWEEN THE NUMBER 
OF UNDETERMINED FUNCTIONS AND THE ORDER OF 
THE EQUATION (iN PARTIAL DIFFERENTIAL EQUA- 
TIONS). 

64. In solving a simple differential equation of the 
first order, we usually arrive at a solution containing one 
undetermined constant. But it is not always so. There is 
a remarkable class of solutions called '' particular solutions/' 
which not only have not, but which cannot have, an un- 
determined constant. Thus, if we express in algebraical 
language the problem " To find the equation to the curve 
whose tangents possess this property, that a perpendicular 
from a given point upon the tangent has a given value ;" 
we have, for general solution, a straight line with one 
undetermined constant; and for particular solution, a 
certain circle, which from the nature of the case does not 

. admit an undetermined constant. The connexion, there- 
fore, between the order of the equation and the number 
of undetermined constants in its solution is not invariable. 

65. But if we assume a formula as solution to a dif- 
ferential equation which we propose to form, the said 
formula containing symbolical constants, we can pre- 
determine that any one of those constants (if the equation 
be of the first order), or any two of those constants (if the 
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equation be of the second order), shall be eliminated from 
the dififerential equation to be formed, and shall therefore 
be undetermined when that equation is solved. Thus, to 
take a very simple instance, let y = aa? +i be assumed to 
be the solution of a dififerential equation, with the constant 

a undetermined. Dififerentiating it, -^" = a; substituting 

this in the given equation, y -x -^ —b = 0, the dif- 
ferential equation required. This process evidently can be 
extended to equations of any order. 

66. It appears thus that the existence of a soluble 
dififerential equation does not necessarily imply an unde- 
termined constant in the solution ; but the existence or 
assumption of an undetermined constant in a solution can 
always be represented by proper form of a dififerential 
equation : and that this theorem applies to equations of 
every order. 

ft 

67. If now we apply similar considerations to Partial 
Dififerential Equations, we are led to the following con- 
clusions. 

68. First, on examining the solution in Article 18, 
&c., it will be seen that, in the cases in which the solution 
of a Partial Dififerential Equation is really efifected, it has 
been done by reducing it to an ordinary integration, just 
as would be done for a Simple Dififerential Equation. This 
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seems to indicate a probability that there may be cases in 
which the solution of a Partial Differential Equation has 
limitations analogous to those of a Simple Differential 
Equation, and therefore has not, and cannot have, an 
undetermined function in its solution. Upon this, how- 
ever, there is not at present any distinct evidence, and 
the matter is indicated as one which appears to merit 
examination. 

69. Second, if we assume a formula containing one 
function <jb(w), where w is a, definite function of x and y 
and where ^ is intended to be undetermined, we can always 
produce a Partial Differeptial Equation of the first order, 
of which the assumed formula will be the solution. 

For we have 

Zy which contains <f>{w) ; 

and we can form 

dz . , 

-J- , which will contain ^(w) and <^'(w) ; 

da . . 

-T- , which will contain ^(w) and ^iyS) ; 

and from these three equations we can eliminate <^(tr) and 

dz 
^{yS)y leaving one equation between z^ ^> y> ^> aJid 

•J- , of the form required. 

D. £. F 
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70i But, third) if we assuble a formula oontaamng two 
faikctioiig ^{w) and '^^)> where ib and 8 are definite foso- 
tions of X and y, and where ^ and y^ are intended to be 
imdetennined, we eannot always produce a Partial Dif** 
ferential Equation of the second order^ of whi<^h the as- 
sumed formula will be the solution. For, we have 

^, which contaiiis <f>{w) and ^(i) ; 

and we can form 

-r- , which will contain <f>{w), 4^'(w), '^(«), y^'{$) ; 
%- , which will eotottfdii ^^\ t^'(^)> ^*)> ^'W ; 

if 

-T-^ , which will contain ^{w\ <l>{w), <l>\to)y '^{s\ 

cPz 
r- , which will contain the same six functions ; 



dxdy 

-7-^ , which will contain the same six functions. 
df 

Here we have dix equations, and no more, with six 
functions on the other side of the equations. We cannot 
therefore eliminate these six functions, so as to leave an 
equation betwe^i 

dz dz 6^ d^z d^z 
^'"^'^'^^ dy' d^' l^y* d^' 
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And thus we cannot always produce the Partial Dif- 
ferential Equation required. The same remark, it will be 
found, applies more strongly to equations of a higher 
order. 

71. It seems not impossible that the failure of at- 
tempts to solve the equations in Article 44 may have some 
connexion with this inability to establish a perfect con- 
nexion between the order of the equation and the number 
of undetermined functions. The subject appears to be 
worthy of greater attention than it has received. 
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